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I describe a procedure for calculating thresholds for quantum computation as a function of er- 
ror model given the availability of ancillae prepared in logical states with independent, identically 
distributed errors. The thresholds are determined via a simple counting argument performed on 
a single qubit of an infinitely large CSS code. I give concrete examples of thresholds thus achiev- 
able for both Steane and Knill style fault-tolerant implementations and investigate their relation to 
threshold estimates in the literature. 
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I. INTRODUCTION 

The threshold for quantum computation is defined as 
the error rate below which, given a number of qubits that 
scales polynomially in the length of the computation, the 
problem size, and the desired fidelity, it is possible to im- 
plement an arbitrary quantum algorithm. Knowledge of 
thresholds is clearly an important tool for use in the de- 
sign of quantum computing architectures, but as of yet no 
simple, unified scheme exists for determining them. The 
great variance of threshold values in the literature is due 
partly to advances and improvements in the field of fault 
tolerance, i.e. error-conscious design, but also partly to 
differing approximations and assumptions regarding er- 
ror models and resources. This makes direct comparison 
of fault-tolerant procedures difficult and tends to limit 
error models to some form of the depolarizing channel. 
The work presented herein was motivated by the desire 
for a method of threshold estimation that could easily 
be applied to a variety of error models and fault-tolerant 
methods. 

One of the keys to achieving fault tolerance is the con- 
struction of low-error encoded states as aids to computa- 
tion. These ancillary states, henceforth termed ancillae, 
provide a baseline against which to check for errors and 
a resource for performing gates that cannot otherwise be 
implemented in a robust way. For the sake of fault toler- 
ance, ancillae are typically prepared in a manner that re- 
duces the frequency of correlated errors at the expense of 
increased overhead in terms of qubits and applied gates. 
Accounting for the vast array of possible methods of an- 
cilla construction has proven a major impediment to the 
design of a generic scheme for threshold estimation. 

Recent work by Reichardt [1] has indicated that, by 
discarding states which show symptoms of error, moder- 
ately sized ancillae can be prepared so well that their con- 
tribution to the failure probability of an encoded circuit 
becomes small. Moreover, my own investigations sug- 
gest that the residual error on ancillae that have passed 
this sort of inspection is primarily due to the verification 
step. When the verification circuit is fault tolerant, this 
implies that the residual error in the inspection basis is 
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uncorrelated. Motivated by these indications, I sought a 
way of determining the threshold given the availability of 
ancillae with independent, identically distributed errors, 
expecting that the result would not differ too much from 
threshold estimates obtained using liberal qubit expen- 
diture during ancilla production. 

This paper describes such a means of threshold calcu- 
lation for fault-tolerant methods employing CSS codes. 
Much of the material presented here can be regarded as 
an elaboration of work by Knill, particularly Refs. [2, 3]. 
As in Ref . [3] , I assume that it is possible to prepare an- 
cillae in logical (encoded) basis states such that errors 
on the component qubits are independent and all qubits 
have the same error spectrum. Using this assumption 
as well as the structure of CSS encoded gates and a few 
error propagation tricks, it is possible to express the fail- 
ure probability of an encoded operation in terms of the 
error probabilities of a single strand of the code blocks. 
In the limit that the number of encoding qubits goes to 
infinity, the criterion for encoded failure becomes partic- 
ularly simple, consequently yielding a threshold in which 
various error probabilities are free parameters, including 
the probabilities of different kinds of Pauli errors on a 
particular gate. 

It should be stressed that the method described here 
does not constitute a constructive procedure for fault- 
tolerant quantum computing. Several details regarding 
syndrome decoding and the selection of the quantum 
code are omitted, and, more importantly, no practical 
prescription is given for preparing the appropriate ancil- 
lae. These caveats do not invalidate the resulting thresh- 
olds, but, presently, the difficulty of obtaining the req- 
uisite ancillae limits their applicability. The algorithm I 
describe is primarily useful cither as a metric for compar- 
ing the performance of fault-tolerant procedures under 
diverse conditions or as a replacement for the numeri- 
cal simulations typically used to generate threshold esti- 
mates. This second application is a particularly apt use 
for my threshold results because the simulations become 
increasingly difficult as the size of the quantum code and 
the proclivity to discard ancillae increase. By contrast, 
the approximations employed in deriving my algorithm 
increase in accuracy with the same parameters. 

The structure of this paper is as follows. Section II 
introduces notation and gives a brief exposition of CSS 
codes, fault tolerance, and thresholds. Section III enu- 
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mcrates my assumptions. Section IV outlines the theory 
needed to generate thresholds using an infinitely large 
CSS code and logical ancillac with independent, identi- 
cally distributed errors. Section V describes how single 
line error probabilities can be determined. In §VI de- 
tails regarding the implementation of my algorithm are 
considered along with the error models to which it ap- 
plies. Section VII performs this analysis for a few cases 
of special interest, including Knill and Steane style fault- 
tolerant quantum computation schemes paired with a se- 
lection of error models. Section VIII presents an alter- 
native algorithm for finite codes. Finally, §IX concludes 
with a brief review of my results and a discussion of pos- 
sible directions for future work. 



II. BACKGROUND 

A. Notation 

Throughout this paper X, Y, and Z are used to denote 
the Pauli operators while X denotes the controlled-NOT 
or XOR gate. I also make use of Hadamard and phase 
gates, which are given in the standard basis by 
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respectively. Together these gates span the Clifford 
group, an important subset of quantum operations. The 
addition of the 7r/4 rotation, 
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completes a universal set for quantum computation. 

I use quantum circuit notation rather extensively. A 
review of this formalism can be found in Nielsen and 
Chuang's book [4]. I deviate slightly from their standard 
in that no decoration is added to distinguish encoded 
circuits from unencoded circuits. 

A group of qubits combined via a quantum code into a 
logical (encoded) state is referred to as a block. Encoded 
states are identified by a line over the contents of the ket; 
encoded gates are similarly decorated when not part of a 
circuit diagram. 

Operations which do not couple qubits that reside in 
the same block are called transversal. This designation 
is in contrast to some of the literature, which insists ad- 
ditionally that transversal operations should apply the 
same component operation to each qubit of a block. 
In this paper, operations satisfying both properties are 
called homogeneous. 

I refer to a single qubit of a block and all qubits that 
directly or indirectly couple to it as a strand. Gates be- 
tween pairs of qubits should be thought of as binding 
then into the same strand, as illustrated in Fig. 1. 

The notation [[n, k,d\] represents a quantum code en- 
coding k logical qubits using n unencoded qubits and 
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FIG. 1: A circuit showing the homogeneous application of 
several encoded gates for the seven qubit code. The inset on 
the right shows a single strand which has been extracted from 
the coding blocks. Note that, modulo a label indicating the 
ancilla's starting state, the strand is identical to the encoded 
circuit. 



having a minimum distance of d. A code with minimum 
distance d can correct any error of weight less than or 
equal to t = ■ All codes dealt with in this paper 

have k = 1. 



B. CSS Codes 

A generalized Calderbank-Shor-Steane (CSS) code [5] 
is a quantum code for which X (bit flip) and Z (sign 
flip) errors can be corrected independently. Canonical 
CSS codes have the additional property that they are 
symmetric with respect to interchange of X and Z. This 
symmetry implies that encoded X, Y, Z, H, and X 
gates can be implemented through transversal applica- 
tion of their unencoded equivalents. As a result of this 
useful property, CSS codes pervade the literature on fault 
tolerance, including this paper. 



C. Fault Tolerance 

Fault-tolerant design is an approach to computation 
and error correction that takes into account the fact 
that every component of a quantum computer is likely 
to be unreliable. Gates, ancillae, measurements, and 
memory are all assumed to err with some probability. 
Fault-tolerant constructions seek to minimize the effect of 
these errors by preventing their spread. This is achieved 
through the use of transversal operations, the expendi- 
ture of qubits, repetition, and teleportation. 

The focus on minimizing the spread of errors is crucial. 
To have any hope of recovering corrupted data, whether 
it be classical or quantum, we need some prior knowledge 
of the kinds of errors that are likely to occur. In the ab- 
sence of any other information, it is generally assumed 
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that errors on different parts of a computer are uncorre- 
cted. Fault-tolerant operations are designed so that the 
failure of a single component does not result in an er- 
ror on many different parts of an encoded state, thereby 
avoiding the creation of correlated errors. 



D. Thresholds 

No matter how skillfully constructed the fault-tolerant 
procedure, there remains for any finite code a nonzero 
probability that too many independent errors will oc- 
cur in a computational step and our data will become 
irreparably corrupted. Consequently, the probability of 
failure approaches one as the length of the computation 
increases. To ensure that we can perform a computation 
of arbitrary length we need a way of making the proba- 
bility of an uncorrectable set of unencoded errors, i.e. an 
encoded error, arbitrarily small. 

The most common method of achieving an arbitrar- 
ily low encoded error rate is known as concatenated 
coding. In concatenated coding, the process of encod- 
ing is divided into many levels. Zeroth-level (physical) 
qubits are used to encode first-level qubits, these first- 
level qubits are then in turn used to encode second-level 
qubits, second-level qubits are used to encode third- level 
qubits, and so on. The basic idea is the following: "If 
encoding qubits reduces the effective error rate then en- 
coding the encoded qubits should reduce the error rate 
even more." This provides us with a plausible sounding 
way of achieving an arbitrarily small error rate; we simply 
add layers of encoding until our error rate is acceptable. 

But encoding will not always decrease our error rate. 
It is possible for our hardware to be so error prone that 
the process of applying an encoded gate and error cor- 
recting is less likely to succeed than simply applying the 
unencoded gate. From this observation arises the idea of 
a threshold error probability pth for quantum computa- 
tion. The threshold is the unencoded error probability 
below which we can achieve an arbitrarily low encoded 
error probability using a number of qubits that scales 
polynomially in the size of the problem. Put another 
way, it is the error probability below which we can com- 
pute indefinitely 

Determining the threshold exactly for a given set of 
assumptions has proven to be a hard problem, but we 
can get some idea of its value through bounds and es- 
timates. Upper bounds resulting from proofs of classi- 
cal simulability have pushed the depolarizing threshold 
below 50% [6, 7], while lower bounds on the order of 
10~ 5 [8, 9] have been rigorously proven through concate- 
nation of explicit fault-tolerant constructions. This paper 
focuses on estimates of the threshold, which have recently 
settled near the 1% mark [1, 2]. 

Estimates of the threshold for quantum computa- 
tion are generally made by analyzing a particular fault- 
tolerant implementation using a specific finite code un- 
der concatenation. Fundamentally, these estimates de- 
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Intuitively, this makes sense; wc would not expect er- 
ror correction to be advantageous when an encoded gate 
or qubit is more likely to fail than an unencoded one. 
Nonetheless, careful consideration of this justification re- 
veals some difficulties with the argument. There are 
many sorts of errors, and it might well be the case that 
the encoded error rate increases for some of them, but 
not all. Even were there only one kind of error, however, 
the equation above would not provide us with a lower 
bound on the threshold, but only an estimate. While, in 
that case, an increase in the error rate at the first level of 
concatenation implies that subsequent layers of concate- 
nation also increase the error rate, the converse is not 
true. To see why, assume we have some fault-tolerant 
procedure for which the encoded failure rate is less than 
the unencoded failure rate. At the first level our code is 
constructed of unencoded qubits that are either perfect 
or have failed. At the second level of encoding, however, 
our code is constructed of singly encoded qubits that may 
be perfect, insufficiently corrupted to result in failure, or 
failed, yet only the last case is considered an encoded er- 
ror. In some sense, the qubits that have not failed are 
now of lower quality than they were at the previous level. 
Thus, the fact that encoding worked at the previous level 
does not guarantee that it will work at the current one. 

Details such as these imply that most threshold results 
are estimates rather than rigorous lower bounds. The 
interested reader can find a more thorough discussion of 
this issue in work by Svore and others [10, 11]. 



III. ASSUMPTIONS 

Some large part of the variance in fault-tolerant thresh- 
old calculations is due to the variety of assumptions em- 
ployed by various authors. In an effort to combat con- 
fusion, my assumptions are listed below in roughly the 
order of decreasing novelty. 

1. All ancillary qubits have independent, identical er- 
ror distributions. 

2. There are no memory errors. 

3. X is the only two-qubit gate. 

4. Any pair of qubits can interact via a two-qubit gate. 

5. Error operators are trace preserving and lack sys- 
tematic coherent terms. 

6. Gate failures are uncorrelated. 

7. Classical computation is freely available. 
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Assumptions 1, 4, 5, 6, and 7 are necessary for my analy- 
sis; the others are convenient but optional. Assumption 2 
obviates the need to consider questions of parallelism, 
gate timing, and the speed of classical computation, while 
Assumption 3 reduces the number of cases that must be 
considered. 



IV. THRESHOLDS FOR HOMOGENEOUS 
METHODS 

Two important observations from §11 provide the foun- 
dation for my method of threshold calculation. The first 
is that fault-tolerant procedures for CSS codes are, to a 
large degree, transversal. The second is that, for the kind 
of CSS codes typically employed, these transversal oper- 
ations can be implemented by applying the same gate to 
every qubit in a block. The sum of these observations is 
that most operations performed in a fault-tolerant proce- 
dure consist of doing the same thing to each of the qubits 
in a block. If one could arrange for all operations to have 
this property, which I refer to henceforth as homogeneity, 
analyzing the behavior of fault-tolerant circuits would be 
greatly simplified. 

Three components of the typical fault-tolerant method 
stand in the way of full homogeneity: ancilla produc- 
tion, syndrome extraction, and recovery. Starting from 
the perspective of threshold estimation, this section ad- 
dresses each of these aspects, partly by keeping in mind 
that the eventual goal is to model errors, not compu- 
tation. Ultimately, the method derived is meaningful 
both as a form of threshold estimation and as a threshold 
bound for idealized resources. 



A. Ancillae 

The ancillae used in CSS-code fault-tolerant proce- 
dures are typically prepared in highly entangled states, 
i.e. in logical basis states. By definition, entangled states 
cannot be constructed without the interaction of the con- 
stituent parts, so there is no a priori reason to think that 
the qubits composing an ancilla will have cither indepen- 
dent or identical error distributions. For fault-tolerant 
procedures, however, the production of an entangled an- 
cilla is usually followed by a homogeneous verification 
circuit, and, in my experience [12], most of the residual 
error probability (of the kind tested for) arises during 
this verification step. With this is mind, I approximate 
ancillae as having uniform error distributions. 

It is important to realize that this assumption is far 
less innocuous than it sounds. Implicitly, I am assuming 
that ancillae of the desired size can be constructed for 
use in a verification circuit, but in subsequent sections I 
take the limit n — > oo. Thresholds given in this limit are 
only practically achievable if an efficient procedure ex- 
ists to prepare logical ancillae. To be efficiently scalable, 
however, a construction routine must have nonvanishing 



probability of generating an ancilla that has good fidelity 
with the desired state. Fault-tolerant schemes using con- 
catenated codes provide a method of achieving this for 
arbitrarily large ancillae, but as a side effect of universal 
quantum computation. At present, there is no known 
method of preparing a logical qubit encoded using a CSS 
code of arbitrary size that does not depend on the abil- 
ity to perform universal quantum computation. Thus, 
absent an explicit recipe for ancillae preparation, the al- 
gorithm presented in this paper does not constitute a 
constructive procedure for achieving any threshold. 



B. Error Location 

Current techniques for locating errors require perform- 
ing a complicated and distinctly non-homogeneous func- 
tion on the output of ancilla measurement. But while this 
classical processing requires knowledge of all the mea- 
surements, its effect, assuming that no more than the 
correctable number of errors has occurred, can be de- 
scribed in terms of the individual qubits. So long as the 
total number of errors present on a measured ancilla is 
less than half the minimum distance 1 , the effect of the 
classical processing is to determine a subset of the mea- 
sured bits that can be flipped to yield an undamaged 
codeword. For the purposes of error correction, know- 
ing this string is equivalent to knowing the location of 
all the errors. While the second kind of information is 
not directly available to a quantum computer, it is quite 
accessible to a theorist treating errors probabilistically. 
I can therefore model the effect of classical processing 
in two steps. First, I determine whether too many errors 
have occurred on a block to permit proper decoding, and, 
if this is not the case, I treat the location of bit flips on 
the measured qubits as revealed. 

I have reduced the process of error location to a non- 
homogeneous failure check and an arguably homogeneous 
revelation step. For a Monte-Carlo simulation, the failure 
check would consist of polling all of the other qubits and 
counting up the number of errors that have occurred to 
see whether they exceeded half the minimum distance 
of the code. If instead we performed a probability flow 
analysis, the expected probability of passing the check 
would be simply 

t=o ^ ' 

where n is the number of qubits in the block, t is the 
maximum number of errors that can be corrected with 
certainty by the code, and pc is the probability that a 



1 Some higher weight errors will also be correctable, but I only 
lower my threshold by ignoring them. 
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particular qubit has an X error at some location (step) 
£, here chosen to be just after the time of measurement. 

Equation (4) suggests a way of recovering full homo- 
geneity. Letting r — t/n, in the limit of large n, Eq. (4) 
becomes 

£ P 0z:) = < 1 ., ^ (5) 
y 1 if p £ > r 

which is again homogeneous from the perspective of a 
simulation. As an added benefit, it is no longer neces- 
sary to concatenate many layers of coding to achieve a 
rigorous threshold; instead a vanishing error probability 
is achieved as the limit of a very large code. This alter- 
native to concatenation is known as large block coding 
or, simply, block coding. 

The preceding paragraphs demonstrate that, for ho- 
mogeneous (independent, identically distributed) errors, 
whether or not an encoded state on a large number of 
qubits fails is determined by the error probability of an 
individual qubit. For this result to be useful, an infi- 
nite family of CSS codes with nonvanishing fractional 
minimum distance must exist. Fortunately, it has been 
shown [13, 14] that CSS codes exist such that r w 5.5% 
for asymptotic values of n. It is not known whether a 
similar claim can be made for CSS codes in which the en- 
coded phase gate can be implemented transvcrsally, but 
this convenience is not necessary for my construction. 

The analysis of this section assumes minimum distance 
error correction, but an identical result applies to er- 
ror correction up to the channel capacity. Gottcsman 
and Preskill [15] have shown that families of general CSS 
codes exist that are asymptotically capable of correcting 
errors up to r as 11%. Hamada [16] has shown that this 
result applies to CSS codes with X-Z exchange symme- 
try as well. 

The appropriate choice for r depends on the purpose of 
the calculation. When the goal is to estimate the thresh- 
old that would be obtained by running a Monte-Carlo 
simulation of a minimum distance decoder, r should be 
chosen to be 5.5%. To obtain the largest bound on the 
threshold for homogeneous ancillae or for comparison to 
threshold estimates that use the channel capacity, it is 
best to choose r = 11%. In other cases it may be desir- 
able to choose a value of r specific to a family of quan- 
tum error correcting codes with special properties, such 
as ease of syndrome decoding or the possession of low 
weight stabilizer operators. 

C. Recovery 

Having diagnosed the location of our errors, the obvi- 
ous way of dealing with them is to apply to each qubit the 
gate which reverses its error. Such a recovery operation 
is inherently inhomogencous since not all qubits will be 
in error, and thus not all qubits will have recovery gates 
applied to them. There arc a number of ways to deal 



with this problem, but I follow the lead of Knill [2, 3] 
and dispense with recovery altogether. I can get away 
with this because a string of Pauli errors can either be 
thought of as an error or as an operator shifting us into 
a different (but equally viable) code space. Put another 
way, we can ignore any errors that we know about since 
we know how to determine the effect (see § V A) they will 
have later in the circuit and this effect is easily accounted 
for when analyzing the results of measurements. 

By a similar argument I need never apply any Pauli 
gates, including those used to implement encoded Pauli 
gates! 

V. ERROR COUNTING 

In the previous section we saw how to modify fault- 
tolerant procedures based on CSS codes so that they 
are fully homogeneous. The advantage of doing this is 
that the error probabilities of qubits within an encoded 
block then become independent and identical. Since each 
strand is functionally identical, it suffices to determine 
the error spectrum for one of them; the probability of an 
encoded failure at any point can be predicted from the er- 
ror probabilities of an individual strand of the transversal 
procedure. As the number of encoding qubits becomes 
large, the fraction of qubits with a particular error ap- 
proaches the expectation for that error. In the limit that 
n — > oo, we can say for certain whether our procedure 
fails on any given step since, in that limit, the probabil- 
ity of an encoded failure becomes a step function. Thus, 
the threshold is completely determined by the probability 
of an encoded failure, and the probability of an encoded 
failure is completely determined by the error probabil- 
ity of a single strand of the blocks. Therefore, in order 
to calculate the threshold I need only determine the er- 
ror probability on a single strand at every point in the 
fault-tolerant circuit. This can be accomplished through 
a combination of error propagation and exhaustive book- 
keeping which I describe in the following subsections. 

A. Pauli Error Propagation 

Pauli error propagation relies on the fact that the Pauli 
group is invariant under conjugation by Clifford gates. 
This implies that any string of Pauli gates followed by 
a Clifford gate is equivalent to the same Clifford gate 
followed by some (possibly different) string of Pauli gates. 
Consequently, it is possible to shuffle Pauli errors to the 
end of a Clifford circuit, thus yielding a perfect outcome 
modified by the resultant Pauli operators (see Fig. 2). 

Furthermore, as explained in appendix A, fault- 
tolerant procedures based on CSS codes need never apply 
non-Clifford gates to the data qubits. Encoded gates out- 
side of the Clifford group are performed by teleporting 
the data onto a specially prepared ancillae. Since only 
Clifford gates are applied, error propagation can be used 
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FIG. 2: Examples of Pauli propagation. Equalities are up to 
an overall phase. 



to determine the effect of a given Pauli error at any later 
point in the procedure. 



VI. PRACTICALITIES 

While the previous section presented the basic algo- 
rithm for determining whether an encoded failure occurs, 
this section deals with details of the error model and the 
implementation that must be considered in any actual 
application of the method. Section VIA identifies the 
error models which can be adapted approximately to my 
method, a question which has been considered for fault 
tolerance in general by several authors [17, 18, 19]. Its 
conclusion, drawn by Preskill in reference [19], is that, 
roughly speaking, coherent errors with random phases 
add like stochastic errors. The finer points of how I im- 
plement the algorithm are covered in § VI B . 



More General Errors 



B. Error Bookkeeping 

Given a gate, say the Hadamard, and a set of proba- 
bilities describing the likelihood of various Pauli errors, 
say p x , p Y , and p z for the errors X , Y , and Z, the post 
gate state can be written as a probabilistically selected 
pure state, such that 



Stochastic Pauli errors are far from the only kind of 
error that can affect a system, but they are an acceptable 
substitute for a variety of other error channels. To see 
why, begin by considering an arbitrary trace-preserving 
error operator £. The action, on a state p, of any such 
error operator can be written as 

£{p) = Y,EiPE] where ^£^=1. (8) 



XH |\&) with probability p* 

YH 1^) with probability p Y 

ZH |\&) with probability p z 

H \*5?) otherwise. 



(6) 



The effect of applying further Clifford gates is to change, 
via error propagation, which Pauli error corresponds to 
each probability, and then to add a second layer of prob- 
abilistic errors. If, for example, we were to apply another 
Hadamard gate our state would become 



X \*S?) with probability p z ■ 

Y |*) with probability 2p Y 

Z \ ^) with probability p x 

|\&) otherwise. 



Px + PyPz 

+ 2p x p z 

- Pz + PyPx 



(7) 

By repeated application of this process it is possible to 
determine the probability of various kinds of errors at 
any point in a circuit composed of Clifford gates. Armed 
with this knowledge we can determine 2 the likelihood of 
an encoded failure or, in the infinite limit, whether an 
encoded failure will happen or not. 



2 For non-transversal circuits the number of terms in our book- 
keeping rapidly becomes unmanageable as the size of the code 



By interspersing errors of this form with perfect quantum 
gates it is possible to model any faulty quantum circuit 
that does not suffer from leakage. When the error oper- 
ators are local, it makes sense to approximate them by 
stochastic Pauli channels. Given a local error operator 
satisfying Eq. (8) , I define the associated stochastic Pauli 
channel to have error probabilities 



Px =]T|tr(^A)| 2 , 

3 

p Y = y^ m e j y )\ 2 , and 

3 

p z =Y J ME ] Z)\ 2 . 



(9) 



increases. 



By design, this channel correctly reproduces the proba- 
bility of measuring that a given Pauli error occurred after 
a single application of the general error operator. Its suit- 
ability in more varied circumstances is the subject of the 
remainder of this subsection. 

In practice, many gates arc required (and therefore 
many error operators act) between each error correction, 
so it is important to know how errors accumulate. As 
for the case of stochastic errors, general trace-preserving 
errors can be separated from the associated circuit pro- 
viding that it is composed exclusively of Clifford gates. 
The separation is accomplished by applying Pauli prop- 
agation to each term in the Pauli-basis decomposition 
of the elements, e.g. Ej, of the error operator. If the 
Clifford circuit is fault tolerant, then error propagation 
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maps single-qubit errors to single-qubit errors (on a given 
encoded block) . Each of the resultant error operators dif- 
fers from the actual error by a simple relabeling of the 
local Pauli basis. Thus, the circuit can be disregarded; it 
is sufficient to consider how the transformed local errors 
accumulate. 

A sequence of s single-qubit trace-preserving errors 
acting on a state p can be written as 



Pex is 



£ s o ■ ■ ■ o £ 2 o £i(p) 

js 32 jl \fe=S / \fc=l / 



(10) 



The probability of measuring, for example, an X error 
on the resulting state is 



E-EE 



32 31 



tr [] E ihk X 



(11) 



By contrast, replacing the error operators with their as- 
sociated stochastic Pauli errors, as defined in Eq. (9), 
yields 



fa = ^^|tr(% fc X)| 2 + 0(p 2 

k=l jk 



(12) 



When p EX and p sx (and the equivalent probabilities for 
Y and Z) agree, the associated stochastic Pauli channel 
is a good substitute for the actual error channel. To the 
lowest nontrivial order in p, the condition for equality 
can be derived as follows. 

Consider each error operator £ as a function of the 
total single application error probability p. Taylor ex- 
panding the elements of £ in ^fp yields 



Eo = I + VPiaoX + P Y + j Z + S I) + 0(p) 
E^o = ^p(ajX + [3 3 Y + l3 Z + 5,1) + 0(p) 



(13) 



where the freedom in the Ej has been used to assure that 
Eq contains the only term independent of p. 

Inserting expanded error operators of the form given 
in Eq. (13) into Eq. (11) and discarding terms of order 
greater than p yields 



££|tr(£; ifefe X)| 5 

k=X j k ^0 



tr 



J E ok X 



\k—s 



(14) 



^tr(^ofc^) ■ 
fe=i 

Thus, to first order in p, the difference between p sx and 



Pex Ps> 



k=l 



Y,^(E ok X) -J2ME ok X)\ 2 

fc=i 

pEi^i 2 ( 15 ) 



E a ° fc 

k=l 



k=l 



=-pEE a ° fca o'- 

fc=l l^k 

As suggested by Preskill [19], this expression has a simple 
interpretation in terms of a 2-D walk composed of s steps 
of sizes ly'paofel- Equation (15) is equal to the difference 
between the square of the displacement for such a walk 
and the expectation of the square of the displacement 
assuming that the walk is random, that is, that stepping 
forward and backward are equiprobable. Thus, the ex- 
pectation of Eq. (15) vanishes if the sign of «ofe is random. 
Taking a slightly different approach, we can treat the en- 
tire expression as the displacement of a 2-D random walk 
composed of s(s — 1) steps of sizes \patok&otyk\- The ex- 
pectation is again seen to vanish when the sign of aok is 
random, but now it becomes clear that the standard de- 
viation will scale like s. Since p EX is also proportional to 
s, this implies that the associated stochastic Pauli chan- 
nel is only really a good substitute when the number of 
qubits being considered is large. 

An identical argument holds for Y and Z errors, show- 
ing that, for my purposes, error models for which the sign 
of the coherent error is random are well approximated by 
their associated stochastic Pauli channel. Conveniently, 
this restriction is preserved under any local relabeling of 
the Pauli bases and therefore applies equally well to the 
original error operators. Examples satisfying the restric- 
tion include all stochastic errors, which have no coherent 
component, and unitary rotation errors where under and 
over rotation are equally likely. Systematic errors, such 
as amplitude damping or a bias towards over rotation, 
are not well modeled, though, in practice, the local rela- 
beling of the Pauli bases imposed by gates will randomize 
these errors somewhat. 



B. Implementation 

In the examples that follow, single-strand error rates 
were determined for three fault-tolerant procedures. For 
each procedure, error rates were calculated for a universal 
encoded gate set, H, X, P, and T gates, as well as for 
an idle step that accounted for the possibility of chang- 
ing the order of X and Z error correction. No checks 
were made on the encoded T gate following its first er- 
ror correction, since the remainder of the gate consists of 
applying the encoded P gate. The encoded P gate was 
assumed pessimistically to be implemented via a telcpor- 
tation process akin to that used for T. For each encoded 
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gate, the maximum was taken over the strand error prob- 
abilities at all measurement steps since in the limit that 
n — > oo encoded failures are caused exclusively by the 
largest relevant error probability at a measurement loca- 
tion. 

The error probabilities for unencoded gates, measure- 
ments, and ancillae were left as free parameters. p r , p AS , 
Pm, Pat, and p BT are used to denote one-qubit, two-qubit, 
measurement, A- type ancilla, and B-type ancilla error 
probabilities where T ranges over the single-qubit Pauli 
errors and AS ranges over the two-qubit Pauli errors. 
Note that ancillae are labeled irrespective of what they 
encode. A-type ancillae are used in locations where Z er- 
rors are more disruptive than X errors, and contrariwise 
for B-type ancillae. In the absence of better informa- 
tion, I assume that A-type ancillae are tested using a 
homogeneous coupling, with discard on failure, for first 
X and then Z errors; the opposite order is used for B- 
type ancillae. In this case I approximate the ancilla error 
distributions as 



Pax 


= Pxz 


rpxi - 


1- Pix + Pxx 


Pay 


= PlY - 


V- PxY 




Paz 


= Piz - 


~Pxz 


(16) 








Pbx 


= Pxi - 


VPxz 




Pby 


= Pyi - 


VPyz 




Pbz 


= Piz - 


- Pxz - 


-Pzi +Pzz, 


first 


order) what 


one would expect if the 



only errors on a verified ancilla were due to undetectable 
errors on the C X gates used to check it. 

It should be emphasized that the ancilla error proba- 
bilities given by Eq. (16) are not the only possible choice. 
They were chosen as a good approximation to the resid- 
ual error following a verification procedure that discards 
the state whenever a problem is indicated. Depending on 
the purpose of the calculation, it will sometimes be more 
appropriate to assign, for example, higher error probabil- 
ities associated with less resource intensive verification or 
different probabilities for different kinds of ancillae. 

The Mathematica program that I use to calculate en- 
coded error rates retains terms up to second order in 
the base error probabilities, but the results given in the 
following sections include only first-order terms. Second- 
order terms were found to be negligible for any plausible 
choice of error model. To understand why, consider a 
simplified error model in which gates can fail in only a 
single way. Let p r be the probability of an individual 
gate failing, and let g\ be the number of gates on which 
a single failure results in an error at location C. Further, 
let <?2 be the number of gates that might participate in 
some pair of failures to yield an error at location C The 



Error 
Model 


Nonzero Error Probabilities 


#1 


n — p n — p n — p 
Pr — 4, Pah — jf- , Pm — 21 

Pax = Pbz = f , Pay = Paz = Pbx = Pby = § 


#2 


Pr = |f , Pah = fg, Pm = 4p, 
Pax — Pbz — y§ , Pay = Paz = Pbx = Pby = ^rf 


#3 


Pah = ^, Pax = Pbz = |§, 
Pay = Paz = Pbx = Pby = §§ 


#4 


Pix = Pxi = Piz = Pzi — f , 
Pax = Pbz = § , Paz = Pbx = f 



TABLE I: Four reduced error models considered in the text 
and in Table III. p r , Pah, Pm, Pat, and p B r represent various 
one-qubit, two-qubit, measurement, A type ancilla, and B 
type ancilla error probabilities where L £ {X, Y, Z} and AH £ 
{I,X,Y,Z}® 2 /{II}. Unspecified probabilities arc zero. 



expected error at location C is then bounded by 



A 9 l h 



92-1 



=9Wr ~ 5l(.92 - l)pi + y (32 - l)pt + 0(pi). 

The inequality arises from the fact that not all pairs of 
failures will necessarily produce an error at C. Even 
ignoring that, however, the second-order terms will be 
negative unless Q2 > 2gi; negative terms may safely be 
neglected since their omission only lowers the threshold. 
Among the examples of the following section, the double- 
coupling Steane procedure, when applied to error model 
#1, has relatively large second-order terms. Yet the 
worst location in that procedure corresponds, roughly, 
to a single-error situation where g± = 13, 52 = 27, and 
p r < .018, for which the ratio of second to first-order 
terms is less than .02. Again, this does not even take 
into consideration the fact that many second-order er- 
rors will be harmless. 



(17) 



VII. SPECIAL CASES 

Having described the operation of my algorithm for 
calculating thresholds, I now apply it to three cases of 
interest. Two of these are variants on a fault-tolerant 
method suggested by Steane [20], while the third case 
is a fault-tolerant telecorrection procedure of the type 
proposed by Knill [2]. 

The unrefined output of this endeavor is the set of max- 
imum strand error probabilities listed in Table II. This 
table specifies a kind of high-dimensional threshold sur- 
face in the space of generic stochastic error models. An 
error channel is below the threshold for a particular pro- 
cedure whenever the maximal strand error probabilities 





Gate 


Maximal Single-strand Error Probability 


Single-coupling Steane 


None 


2p AY + 2p B Y + Piy + 2pM + Pxy + pxz + Pyi + 3p Y x + 3p Y v + Pyz + 4p Z x + 3p Z Y + vaecx(2p AZ + 2p B z + piz + 2p x + 2p Y + 3pyi + Pyx + 2p Y z + 4p ZJ + 

3p Z z , 2p AY + 2p BY + Piy + 2p M + p X Y + Pxz + Pyi + 3p Y x + 3p Y Y + Pyz + 4p Z x + 3p Z Y ) 


H 


p AX + 2p AY + p AZ + Pby + Pbz + Pix + Piy + 2p M + p x + Pxz + 2p Y + 2p Y i + 2p Y x + 2p YY + Pyz + 2p Z i + 3p Z x + 2p Z Y + Pzz + max(p X i + Pyi + 2p Y z + 

pz + Pzy + pzz , Piy + piz + px + pxx + 2p X Y + Pyx ) 


C X 


Pay + Pby + 2p [Y + 3p M + 2p XY + 2p xz + 2p Y i + 3p Y x + 3p YY + 2p YZ + 5pzx + 3p ZY + max(p AZ + 2p BY + 3p B z + 2p rz + 3p x + 3p Y + 3p Y i + 2p YX + 

Pyz + 5p Z i,3p AX + 2p AZ + Pbx + 5pix + 3p rY + 2p X i + 3p X x + Pxy + 2p ZY ) 


T P 


Bhv -+r)j\* -\-Vvt +Dvy -\-T)vy -\-Vv7 -\-T>7x -\-Vvy -+max( 04 x Av -\-T)nx -\-2Vix -\-2Vjy j t'Dm-\-T)xt-\-T)xx -\-T)xv -f A-TI7Y T)by -+2DR7 -\-T)x -+Z>v +077+077 1 


Double-coupling Steane 


NOI16 


Pay + Pby + 3p n - + Pm + 3p X Y + 3p X z + 3p Y i + §Pyx + 7p Y Y + 3p Y z + 5p Z x + §Pzy + max(3p A Y + 4p AZ + p BZ + 3p TZ + Px + Py + 2p Yr + ^Pyz + 5pzi + 

2p Z Y + 7p Z z,p A x + 4p B x + 3p BY + 5pix + 2piY + 3p X i + 7p X x + ^Pxy + 2p Y x) 


H 


2p AY + 2p AZ + Piy + Piz + Pm + Px + Pxi + Pxx + 2p K Y + 3p X z + Py + 3p Y i + ^Pyx + 5p Y Y + §Pyz + 2p Z i + 3p Z x + 4p Z Y + 3p Z z + max(p A x + p A v + 

pix + 2pn' + Piz + Pxx + 2p X Y , Pby + Pbz + Pxi + Py + 2p Y i + 2p Y z + Pz + Pzi + Pzz ) 


C X 


p A Y + Pby + 4:p IY + Pm + 5p X Y + 5p X z + ^-Pyi + 4p Y x + 7pYY + 5p Y z + 4pzx + 4p Z Y + max(p^y + 2p AZ + Pbz + Piy + 5piz + Px + Py + 2p Y z + 4p Z i + 

3p Z Y + 7pzz , p A x + 2p B x + Pby + 4p IX + 5p X i + 7pxx + 2p XY + Pyi + 3p Y x ) 


T, P 


2p BY + Piy + Pm + Pxy + Pyx + 2p YY + Pzx + 2p ZY + max(2p BX + 2p IX + Piy + 2p X i + 3p X x + 2p XY + 2p X z + 2p Y i + 2p Y x + Pyy + 2p Y z + Pzx , 2p BZ + 
Piz + Px + Pxz + Py + Pyi + 2p Y z + Pzi + 2p Z z , p A x + p ay + 3p B x + Pby + 3pix + 2p [Y + 3p X x + 2p XY + 2p Y x + Pyy + 2p Z x + Pzy ) 


Knill 


None, 
T, P 


p A Y + Pby + Pm + Pyx + Pyy + Pzx + Pzy + max(p AZ + p B z + Px + Py + Pyi + Pyz + Pzi + Pzz , p A x + Pbx + Pix + Piy + Pxx + Pxy ) 


H 


Pbx + 2p B Y + Pbz + Pm + Px + Py + Pyx + Pyy + Pzx + Pzy + min(p /x + Piy + Pxx + Pxy , Py + Pyi + Pyz + Pz + Pzi + Pzz) 


C X 


p A Y + Pby + Pm + Pyx + 2p Y Y + Pzx + Pzy + max (p AY + 2p AZ + p B z + Px + Py + 2p Y i + Pyx + 2p Y z + 2p Z i + pzx + Pzy + 2p Z z , p A x + Pbx + pix + 
Piy + Pxi + 2p X x + 2p XY + Pxz + Pyi + Pyx + Pyz , p A x + 2p B x + Pby + 2p TX + 2p TY + 2p X x + 2p XY + Pyx + Pzx + Pzy , p A z + Pbz + Piy + Piz + Px + 

Pxy + pxz +Py + Pyi + 2p Y z + pzi + Pzy + 2p Z z) 



TABLE II: Maximum error probabilities for a single strand of various encoded gates for a selection of fault-tolerant procedures. The subscripted p's refer to probabilities 
of error of various unencoded operations with X, Y, Z denoting Pauli errors, M measurement errors, and A and B ancilla errors. Thus p x , Piz, pzi, Pm, and p AY 
denote the probabilities of a single-qubit-gate X error, a Z error on the target of a °X, a Z error on the control of a C X, a measurement error, and a Y error on an A 
type ancilla. A procedure is below threshold when the maximum error probability for all gates is less than r, the fraction of errors corrected asymptotically. 
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Thresholds for Homogeneous Ancillae (r) 


' — — -_____J3rror Model 
Procedure — — — ______ 


#1 


#2 


#3 


#4 


Single-coupling Steane 


0.15 


0.06 


0.24 


0.29 


Double-coupling Steane 


0.16 


0.10 


0.18 


0.29 


Knill 


0.35 


0.15 


0.50 


0.67 



b) 



TABLE III: Threshold for ancillae with homogeneous errors 
given in units of r, the correctable error fraction, for the 
three procedures and four error models I consider as exam- 
ples. These thresholds were obtained by substituting the pa- 
rameters given in Table I into Table II and requiring that 
single-strand error probabilities not exceed r. 



for that procedure are lower than the fraction of errors 
that are correctable asymptotically. For the purpose of il- 
lustration, however, it is more useful to consider less com- 
plicated error models. Table I defines four reduced error 
models in terms of the generic stochastic error model. 
Since these reduced error models have only a single free 
parameter, their threshold surfaces are simply numbers. 
Table III lists thresholds for three procedures and four 
reduced error models in terms of r, the asymptotic cor- 
rectable error fraction. 

The following subsections provide supplemental infor- 
mation specific to each procedure, including qualitative 
reviews of the procedures, circuit diagrams for encoded 
gates, and commentary on the thresholds given in Ta- 
ble III. 



A. Steane's Method 

Steanc's method is characterized by the extraction of 
error information through the coupling of encoded data 
to ancillae prepared in simple logical states. Using an 
ancilla prepared in the logical state |0) (|+)) it is possible 
via a single transversal two-qubit operation to extract 
information about the location of Z (X) errors on all 
qubits in a data block. 

Typical instantiations of Steane's method employ mul- 
tiple extractions to guard against errors made during the 
coupling and measurement process. Often [1, 21, 22] the 
number of extractions performed is conditional on their 
output. I deviate from this rule by demanding a fixed 
number of couplings. This a sensible choice for my anal- 
ysis since, up to rearrangement of qubits, the output of 
an extraction becomes deterministic as the number of 
qubits in an encoding approaches infinity. Moreover, if 
p is the probability of an error occurring, requiring the 
sequential agreement of j extractions reduces the prob- 
ability of misdiagnosing an error on a particular line to 
order p> . Since I ultimately retain only first-order terms 
I need only consider single and double-coupling Steane 
procedures. 



A A -®-wC 
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Gate 
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H 



X 



P, T 



Circuit 



Fz - Fx - Fx - F z - F z - Fx 



- Fz - Fx - H - Fx - Fz - 



- Fz - Fx 



— Fx — F : 



Fz - Fx - 



Fx — Fz — 



Fx\-\F Z 
A 



Fz - 



FIG. 3: Encoded circuits for the single-coupling Steane pro- 
cedure. Handling of the data is minimized at the cost of 
syndrome verification. Parts a and b display the circuits for 
finding X and Z errors, respectively. Part c lists the circuits 
analyzed to determine the encoded error rates for this exam- 
ple. 



1. Single- coupling Steane Procedure 

Two-qubit gates require a controlled interaction be- 
tween two otherwise isolated quantum systems. Conse- 
quently, they are often the most error prone gates in a 
universal set. In such cases, the factor limiting the prob- 
ability of successful error correction may be the number 
of times that two-qubit gates must be applied to the data 
in order to reliably diagnose errors. For Steane's method, 
this interaction is minimized by coupling to the data once 
per X correction and once per Z correction, as shown in 
Fig. 3. 

Table III shows that, relative to the other procedures 
considered, the single-coupling Steane procedure per- 
forms most strongly for error model #3. This is in line 
with our expectations since model #3 includes only two- 
qubit gate errors and the resultant ancilla errors. Surpris- 
ingly, it also does rather well overall, suffering in compar- 
ison to the double-coupling Steane procedure only for er- 
ror model #2 where measurement errors dominate. The 
single-coupling Steane procedure lacks a means of syn- 
drome verification, so any errors in syndrome measure- 
ment are transferred directly to the data. Nevertheless, 
my results demonstrate that moderate single-qubit and 
measurement error probabilities can be tolerated when 
high quality ancillae are available. 



2. Double- coupling Steane Procedure 

When two-qubit gates are relatively reliable, the dam- 
age done during the extraction of error information can 
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be limited by preparing ancillac such that they include 
few errors capable of propagating to the data. Under 
these circumstances, it is often advantageous to verify 
error diagnoses by coupling to the data more than once, 
as shown in Fig. 4. 

For the double-coupling Steane procedure, error model 
#1 is especially interesting because it was chosen in im- 
itation of the error model used by Reichardt [1] in his 
numerical estimation of the threshold for a Steane style 
procedure on a 49 qubit code. My threshold of 0.90% 
for asymptotic minimum distance decoding is quite close 
to his value of roughly 0.88%. The extraordinary agree- 
ment of these two estimates is a coincidence, as can be 
seen from my discussion of finite codes in §VIII, but their 
rough equivalence illustrates the value of my idealized al- 
gorithm for approximating the encoded error rates used 
in threshold estimation. Implementing Stcane's method 
with 49 qubit ancillac prepared via a liberal discard pol- 
icy yields roughly the same threshold as implementing it 
with very large ancillae prepared such that errors due to 
the two-qubit gates used for verification dominate. 

As expected, relative to the other two procedures, the 
double-coupling Steane procedure performs most favor- 
ably for error model #2. Somewhat surprisingly, how- 
ever, it still undcrpcrforms the Knill procedure. The rea- 
son for this is most easily understood by considering the 
limiting case in which only measurement errors occur. In 
the absence of any other source of error, measurement er- 
rors have no effect on either the double-coupling Steane 
or the Knill procedure until their probability exceeds r; 
beyond that point both procedures fail with certainty. 
Thus, the two procedures cope with measurement errors 
equally well, but the Knill procedure handles other kinds 
of gate errors more effectively. 

Error models #3 and #4 demonstrate small gains in 
the threshold that can result when two-qubit gate errors 
have some underlying structure. Model #3 is a pure two- 
qubit-gate depolarizing error model (which includes the 
associated ancilla errors) while #4 is a model in which 
two-qubit gates malfunction by producing either an X or 
a Z error on either the control or the target. Given the 
highly restricted form of error model #4 it is discouraging 
that the threshold increases by less than a factor of two 
over that of error model #3. Nonetheless, the possibility 
of larger gains in a procedure designed to take advantage 
of some particular error distribution is not ruled out. 



B. Knill's Method 

Knill's method for fault-tolerant quantum computa- 
tion [3] utilizes telecorrection (see Fig. 5). Telccorrec- 
tion, error correction by teleportation, requires the ap- 
plication of even fewer two-qubit gates to the data un- 
dergoing correction than a single-coupling Steane error 
correction. The cost of this innovation is a greater re- 
liance on the ability to produce high quality ancillae. In 
effect, Knill's method exchanges the error distribution of 
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FIG. 4: Encoded circuits for the double-coupling Steane pro- 
cedure. Syndrome information is extracted twice, and qubits 
implicated both times are presumed to be in error. Parts 
a and b display the circuits for finding X and Z errors, re- 
spectively. Part c lists the circuits analyzed to determine the 
encoded error rates for this example. 



an encoded data block for that of one half of an ancilla 
prepared in a logical Bell state. Failure occurs only when 
too many errors are present at the time of measurement 
to correctly identify the encoded Pauli operator needed 
to complete the teleportation. Conveniently, encoded 
singlc-qubit gates can be accomplished by performing the 
teleportation using a logical Bell state prepared with the 
desired gate already applied to one half. 

For my implementation of Knill's method, error model 
#4 achieves the highest threshold, though physical sys- 
tems displaying this sort of error seem unlikely. Error 
model ^1 provides another check of my algorithm, since 
its parameters are also roughly those used by Knill [2] in a 
paper on telecorrection. Setting r to 11% for the channel 
capacity for CSS codes, I find that the threshold for this 
model is 3.9% compared to Knill's estimate of 3% and his 
extrapolation of up to 5%. The approximate agreement 
between these values is satisfying, though an exact match 
is not expected since Knill assumes that errors on up to 
19% of the qubits can be corrected, an assumption that 
derives from bounds on the channel capacity for general 
quantum codes [23]. 

Of course the most striking aspect of Table III is that 
the Knill procedure yields a higher threshold for every er- 
ror model. As with the single-coupling Steane case, this 
derives partly from my assumptions regarding ancillac. 
In particular, Steane's method was designed to utilize 
ancillac for which either correlated X or correlated Z er- 
rors could be minimized, but not both, a situation certain 
to favor his approach. A second but lesser objection can 
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FIG. 5: Encoded circuits for the Knill procedure. Error cor- 
rection is performed by teleporting the data, minus the er- 
rors, using an entangled two-logical-qubit ancilla; the precise 
location of errors is unimportant so long as the result of the 
logical-qubit measurement is correctly decoded. Parts a and 
b display the circuits for correcting errors when the output 
of the previous step was dominated by Z and X errors re- 
spectively. Part c lists the circuits analyzed to determine the 
encoded error rates for this example. Ancilla error distribu- 
tions are used as the input to these circuits since that is the 
only remaining source of error after a successful teleportation 
correction. 



be made that I set the ancilla error probabilities equal 
for all gates and all methods, ignoring the fact that some 
methods, such as Knill's, and some gates, such as the T 
gate, will require more complex ancillae which may in 
turn be more error prone. Substantially more detailed 
ancilla information would be needed to evaluate the im- 
portance of this effect, but the overall character of my 
results is unlikely to change since that would entail in 
excess of a two-fold increase in the error probabilities for 
logical two-qubit ancillae over those for ancillae prepared 
in a singlc-qubit logical state. Thus, so long as resource 
considerations do not limit our ability to discard suspect 
ancillae, and therefore to make very high quality ancillae, 
Knill's method will provide the highest thresholds. 



VIII. FINITE CODES 



and its interpretation as an idealized threshold bound is 
completely lost. Fundamentally, the complications that 
arise are all due to the fact that the success or failure 
of various portions of an encoded gate are no longer de- 
terministic. This section explains how to deal with the 
associated difficulties and concludes with a brief demon- 
stration of the algorithm for a [[49, 1, 9]] code. 

In the examples of §VII, I establish a background er- 
ror rate by performing an initial error correction, but for 
finite n this initialization is not guaranteed to succeed. 
Though the failure of the initial error correction is prop- 
erly assigned to the previous encoded gate, the residual 
errors will differ dramatically depending on whether it 
occurred. This presents no problem when only a single 
level of encoding is employed since any encoded failure is 
considered a failure of the computation. In concatenated 
coding schemes, however, failed encoded qubits are cor- 
rected at higher levels of encoding. Their continued use 
is problematic since an encoded gate failure may be cor- 
related with subsequent encoded failures. Nevertheless, 
I recommend calculating the encoded error rate for finite 
codes using the assumption that the initialization did not 
fail, a choice that requires no modification to the case for 
large n. 

Likewise, calculation of the single line error rate pc 
proceeds without modification. For finite codes, how- 
ever, the maximum tolerable single line error rate be- 
comes a nontrivial function of the encoded error rate that 
we wish to achieve. The probability of an unrecoverable 
error never goes to zero, so it is necessary to perform the 
summation in Eq. (18) to determine the portion of the 
encoded error rate due to any particular location. 

The possibility of failure must be considered at many 
points in the circuit since statistical fluctuations will pro- 
duce unrecoverable errors at a variety of locations. Typ- 
ically, encoded failure probabilities at various locations 
will be strongly correlated, but the exact nature of these 
correlations is difficult to predict. Thus, the best I can 
do is to bound the encoded failure probability, 



max E{(pc) < 



Encoded 
error rate 



< 



E 

ces 



Etipc) (19) 



Prior to taking the limit n — ► oo, the expression for 
the probability of an encoded error at a location £ was 



Ei(pi 



E 

i=t+l 



Pc^-PcY 



(18) 



where t is the number of correctable errors and pc is the 
probability of a relevant error on a single qubit at the lo- 
cation in question. Using this expression, the programme 
of §IV can be implemented for finite n. In doing so, how- 
ever, the simplicity of the algorithm suffers somewhat, 



where S ranges over the locations of every post- 
initialization output, that is, syndrome measurements 
and the final state of the data with regard to both X 
and Z errors. 

To clarify the changes outlined above, consider the ex- 
ample of the double-coupling Steane procedure imple- 
mented using a [[49,1,9]] quantum code and subject to 
the error channel defined by error model #1. For the 
encoded C X gate, the set of single line error probabili- 
ties corresponding to X errors at the eight locations of 
post-initialization syndrome measurement and X and Z 
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errors at the two output locations of the data is 

f 47p 43p 43p Alp 39p 33p 37p 
IPs} - | — , 

31p 9p 9p 3p 3pl 

~8~'T'T'T'Tj ' 

Solving Eq. (19) subject to the restriction that the en- 
coded error rate is exactly p yields solutions in the range 

.0036 > p > .0023. (21) 

Repeating this process for each of the other encoded 
gates and taking the minimum over the upper and lower 
bounds produces a threshold of 

.0034 > Pth > .0023 (22) 

where, of course, the caveats discussed in §11 D regarding 
concatenated threshold estimates all apply. This exam- 
ple provides a particularly apt comparison to Rcichardt's 
threshold estimate for the [[49, 1, 9]] code [1]. The thresh- 
old calculated here is roughly a third of that estimated 
by Reichardt. The difference presumably springs from 
the superiority of his rule for syndrome extraction when 
n = 49. 



IX. CONCLUSION 

I have described a general algorithm for generating 
thresholds provided that ancillae with independent, iden- 
tically distributed errors are available as a resource. My 
approach applies to most fault-tolerant procedures em- 
ploying CSS codes. It relies on the fact that nearly all 
elements of such a procedure are homogeneous, that is, 
transversal with identical components. Inhomogcneous 
elements are either eliminated, as for classical syndrome 
processing and the application of recovery unitaries, or, 
in the case of ancillae, replaced with homogeneous equiv- 
alents. This allows me to calculate the probability of fail- 
ure for encoded gates in terms of the error probabilities 
associated with a single strand of the encoded blocks. In 
the limit that the number of encoding qubits approaches 
infinity, a criterion for success becomes simply that the 
probability of finding an error never exceed the fraction 
of the encoded qubits on which said error can be cor- 
rected. When this is satisfied, it is possible, in the limit 
of infinite block size, to compute indefinitely, and our 
base error rates are, by definition, below threshold. 

The value of considering thresholds for homogeneous 
ancillae is that they can easily be calculated for a variety 
of fault-tolerant procedures and error models, thereby 
providing a relatively simple metric for comparison. Sec- 
tion VII includes thresholds for computation for three 
fault-tolerant procedures and four error models. One 
of the procedures considered is based on a method of 
telecorrection used by Knill, while the other two are vari- 
ations, in that the number of syndrome extractions is 



fixed, on Steane's approach to achieving fault tolerance. 
The error models considered are a full depolarizing error 
model, a depolarizing error model with increased mea- 
surement errors, a depolarizing error model for two-qubit 
gates exclusively, and a restricted two-qubit-gate error 
model. Holding the total probability of an error con- 
stant, small improvements are observed in the threshold 
for certain choices of the two-qubit-gate error model. For 
the procedures examined, the threshold increases by less 
than a factor of two, but larger gains may be achievable 
using fault-tolerant procedures tailored to a specific error 
model. With regard to comparisons between procedures, 
the single-coupling Steane procedure is shown to out- 
perform the double-coupling procedure when two-qubit 
depolarizing errors dominate, but the double-coupling 
Steane procedure does substantially better when mea- 
surement errors are likely. I also find that Knill's ap- 
proach outperforms that of Steane for all error models 
considered, a conclusion that is likely to hold so long 
as correlated ancillary errors are rare and the ancillae 
needed for Knill's method are not appreciably more er- 
ror prone than those employed by Steane. 

Idealized thresholds aside, my algorithm is useful as 
a means of approximately computing the logical error 
rate for a single level of encoding, which is an established 
method of estimating the threshold for quantum compu- 
tation. The two treatments yield similar outcomes be- 
cause numerical estimates of the encoded error rate typi- 
cally prepare ancillae in a way that maximizes their qual- 
ity at the cost of additional resource overhead. Ancillae 
prepared in this manner have error distributions approx- 
imating my ideal of independent, identically distributed 
errors. The basic algorithm uses the infinite limit to ob- 
tain simple analytic results, but an alternative (and less 
rigorous) algorithm for finite codes is described in §VIII. 
Both methods were shown to yield results in rough ac- 
cordance with the depolarizing threshold determined by 
Reichardt [1] for the [[49,1,9]] code. For telecorrection, 
in the limit n — > oo, my estimate of the depolarizing 
threshold was consistent with the range of values deter- 
mined by Knill [2]. 

I provide Table II for researchers who wish to estimate 
the threshold of the Steane or Knill method given a par- 
ticular error model. For those interested in analyzing 
different fault-tolerant methods, my Mathematica pro- 
gram is available at http://info.phys.unm.edu/~beastin 

Much further work remains to be done on this subject. 
One topic of interest is the tailoring of fault-tolerant pro- 
cedures to the error model. It was my initial hope that 
restricted error models would greatly improve the tol- 
erable error rate. While this was not observed in exist- 
ing fault-tolerant procedures, the possibility remains that 
tailored procedures could improve the threshold. A sec- 
ond possibility is the extension of my analysis to include 
memory errors, which promises to be a straightforward, 
if unbeautiful, endeavor. The most valuable addition, 
however, would be to explicitly define methods of ancilla 
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construction and determine the degree to which they dif- 
fer from my ideal. Constructing ancillae to my specifi- 
cations is an extremely difficult problem, but one whose 
solution would have a strong impact on the theory of 
quantum computing in general and this work in particu- 
lar. A scalable method for producing ancillae with inde- 
pendent, identically distributed errors would enable the 
algorithm presented here to be employed for the calcula- 
tion of rigorous lower bounds on the threshold without 
any caveats about idealized resources. 
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APPENDIX A: tt/4 ROTATION 

The 7r/4 rotation, T, is not a member of the Clifford 
group. This means that T does not take Pauli strings to 
Pauli strings under conjugation. A logical T gate would 
also have this property, thereby confounding the simple 
error propagation routine used in this paper. For CSS 
codes, however, the encoded T gate is not applied di- 
rectly; instead, the logical state 
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is prepared and the data qubit is effectively teleported 
under the T gate using the following encoded circuit: 



|e) 



px - 



(A2) 



If necessary, P can be implemented in the same way as T; 
otherwise the gates applied are Clifford operations (both 
in terms of the encoding and its constituent qubits), so 
error propagation proceeds without a hitch. 

There is a modest sleight of hand here in that we would 
not, if we so desired, be able to apply error propagation 
during the construction of |<9). For my purposes this 
is unimportant since I do not model the construction of 
ancillae. 



